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ABSTRACT 

The problem of determining atomic lifetimes from bean- 
foil decay curves which are subject to cascading effects is 
discussed. The ANDC method of analysis, which in theory can 
exactly account for all cascade effects, is described and a 
practical method of applying it is developed. Several 
quantities which result from the analysis, and the possible 
Variations on the method of analysis, can be used as self- 
consistency checks on the applicability of the method toa 
given situation. 

The resuits of a test on three decay situations using 
synthesized data are given along with examples of its 
application to real data. Using these results, several 
rules for applying the ANDC method and interpreting the 
results are developed. Recommendations for obtaining 
experimental decay curves most suited for ANDC analysis are 
also mades The computer program developed to perform the 


analysis is included along with instructions for its use. 
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CHAPTER <I 
INTRODUCTION 
An important complement to any theory is its 
verification by experiment. The application of quantum 
mechanics to many-electron atoms is not a straightforward 
task. Many assumptions and approximations are necessary in 
order to calculate atomic energy levels.! Before the advent 
of beam-foil spect roscopy2, hereafter referred to as BFS, 
the amount of experimental data available for comparison was 
relatively small due to the limitations of the conventional 
sources {Shock tubes,electric discharges and furnaces). The 
Main advantages of BFS are that one can, in principle, 
excite almost any ionic state of any element and that one 
can also follow the decay of that state with time. This 
enables one to study trends in oscillator strengths along 
various seguences (e.g., isoelectronic and homologous). 
Comparisons of these experimentally determined trends with 
theoretical predictions provide a good test of the accuracy 
of the approximations used in the quantum mechanical 
calculations (CCF 76). Accurate oscillator strengths for 
neutcal and highly ionized atoms are also needed for the 


determination of astrophysical quantities (e.g, particle 


1 Theorists work in terms of transition probabilities Ajj of 
oscillator strengths f;;, while BFS experimentalists 
determine lifetimes. The relationships between the three 
quantities are given in section 1.1 below. 

2 Pioneered by S. Bashkin in 1964 (Ba 64). See also the 
recent review article on BFS (Be 77). 
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densities) and for studies ees and plasma physics. 

The past ten years have seen four international 
conferences devoted entirely to BFS.! Major topics at these 
conferences have been the determination of atomic lifetimes, 
and possible methods to improve the accuracy of these 
determinations. It should be mentioned that BFS has also 
been used to determine fine and hyperfine structures, Lande 
g-values and for general spectral analysis of ionized atoms. 
Examples of such studies may be found in these conference 


reports. 


A mono-energetic beam of ions, having velocity 
typically in the range from 1 to 5 mm/ns, is magnetically 
analyzed for the desired mass and directed through a thin 
foil (usually of carbon and having a thickness of typically 
500A). The interaction with the foil results in a beam of 
ions at various levels of excitation and ionization, 
depending on the beam energy. The excited states then decay 
by the emission of radiation and the intensity of a 
particular transition can be measured as a function of time 
(see Fig. 1.1). 

The basic data recorded is a series of triplets 
(dj,Nj,T;), where d; is the distance from the foil and N: is 
the number of photons accumulated in time T;. Time is 


1 Proceedings: Ba 68,MBB 70,Ba 73,SP 76. 
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Fig. 1.1 Basic beam-foil experimental 
configuration (not to scale). 


Spectrometer 


determined by a pre-set level for integrated beam current as 
measured by the Faraday cup, i.e., intensity is normalized 
by beam current. 

In the past, this photon counting procedure! has been 
the most common method of collecting beam-foil data. 
However, some experimental groups, including the one at the 
University of Alberta, are now using multi-channel analysis 
techniques to obtain more reliable data with greater 
efficiency. In the multi-channel analysis(MCA) method 
several sweeps (varying distance from the foil) are summed 


together to obtain the final decay curve. The main 


1 The general experimental design and procedures that have 
been used at the University of Alberta are described in IL 
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advantage of this method is that during a run any 
Significant change in the decay curve of the latest sweep, 
aS compared to previous sweep(s), indicates a change in the 
excitation process (e.g., beam instability, change in foil 
condition) and the run can be terminated at this point 
Without a significant loss of data. Thus the MCA method 
allows for the constant monitoring of the quality of the 
data and minimizes the time and foils spent in acquiring the 
data. 

Obviously a primary condition for determining accurate 
lifetimes is the acquisition of the best decay data 
possible. However, the means of analyzing this data once it 
has been obtained is a separate problem and is the topic of 
this work. Thus details of experimental equipment and 
procedures wiil not be discussed. 

If only one atomic level was excited, one could observe 
the exponential decay in intensity and thus measure the 
lifetime, tT, of that level directly from the slope of the 


logarithmic decay, since 


T 55 (t) a -dN,/dt and 


or In(T, 5) == / tb conse. 


where N. is a constant corresponding to the count which 


0 
would be recorded at the foil position. In the event that 


the transition i+j is the only (or dominant) mode of decay 
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fromelevel 1, then*t;°is simply V/Ri ge where Ais is the 
transition probability. 
Fig. 1.2 Simple two level system. 
5 
j 
This is related to the oscillator strength, f. by the 
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relation 
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where Jivdj are the statistical weights of the two levels. 
If more than one decay mode can significantly depopulate the 


level i, then TAS) A ss Conversion from To tatop ie, then 
J 


ji 
requires a knowledge of the relative values of the Aya 
which can be obtained from branching ratio measurements or 


from theoretical calculations. 


1.2 Problems of the 
The wide range 


of a mixed blessing 


the excitation process. 


Method 
of states available with BFS is somewhat 
because of the non-selective nature of 


Levels above the one under 


consideration are also excited and decay into the level of 
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interest, repopulating it and masking the simple exponential 
behaviour of a single decay. This problem of cascading 
levels is the major difficulty in determining lifetimes. 
Common methods of dealing with the problem are described 
briefly in the next section. 

Occasionally, two different transitions have similar 
wavelengths so that the lines overlap. This problem always 
occurs to some degree, but in most casSeS iS not a 
Significant source of error. 

The beam density is low (~10+5cm-3) so that trapped 
radiation, stimulated emission and ionic recombination 
effects are negligible (these are usual problems with many 
conventional sources). However, this means that the 
intensity is correspondingly low. Also, the foil ages and 
eventually breaks due to the action of the beam, giving a 
finiter tine in which one ‘can collect data. This limits BFS 
to relatively strong transitions unless multi-channel 
analysis techniques are used (PWVL 77). The time scale is 
determined by the velocity of the beam, thus the accuracy of 
the Lifetime is limited by the uncertainty in the velocity. 
Various techniques can be used to obtain the beam velocity 
to within 1%, which is usually less than the statistical 
error resulting from the cascading problem to be discussed 


in the next section. 


In order to take into account the presence of 
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cascading, one must first examine its contribution to the 
decay curve. The simplest case which shows the effect of 


cascading is the three level system shown in Fig. 1.3 . The 


Fig. 1.3.Three level decay scheme. 


a 


intensity of the primary decay is T51(t)=a2N5(t) but, 


dN. (t) a 
CADE = —a2N. (t) 52 asN,(0)e 03 ’ 


dt 


therefore 


(oo = Oia Je 


d(N. (t)e%?*) 
__ Oh asN.(O)e : 


dt 


Integrating both sides of the equation with respect to time 


gives: 


Cintas OL 3 (a2-03)t 


QA27A3 


N,(0)e + const. 


Setting t=0, one finds that 
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const = N, (0) ayer N3(0) . 


Substitution of this gives the population of the primary 


level: 


- O 3 =O5.C O 3 -a3t 
N, (t) | Np (0) + =o N3 (0) | e + —2_N,(0)e 


Q 3 Obey hb) 


Thus the primary intensity is given by the expression: 


To, {t) 


2 


O 3 =O 9 t 
EAE N,(0)e . 


3 Q3 —-aot 
| Np (0) + Nn, (0)| s to 


O32 -6 
This shows a particular instance of the general problem: the 
decay curve has as many exponential terms as levels which 
cascade into it plus the one due to itself.! Usually a few 
of these terms dominate; however, this still leaves some 
difficulty in obtaining the desired lifetime. 

The most obvious thing to do is to try to fit the decay 
curve to a function consisting of a sum of exponential 
terms. Usually it is not worthwhile or necessary to try 
more than three exponentials plus a constant background. 
This has its problems, however, since the fit is not unique, 
and one must make a careful analysis to find the "correct" 


parameters. Computer programs have been written to do this 


1 The general resuit for an arbitrary decay scheme is given 
in Appendix I. 
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analysis (e.g., HOMER see IL 74) and has been found to be 
adequate in most cases. The estimated errors in the 
parameters obtained from such fitting procedures are only 
valid if the form of the fitting function is correct (PDKI 
78), and determining the correct form can be difficult when 
severe cascading is present. 

Integral and differential techniques are useful when 
the primary lifetime is much shorter or longer than any of 
the significant cascades. If one differentiates the 
intensity curve, the relative amplitude of the term with the 
shortest lifetime is increased. Similarly, if the intensity 
curve is integrated the longest lifetime term increases in 
relative amplitude. This may be seen by the following 
example of an intensity curve with two exponential terms 


corresponding to lifetimes of two and ten seconds. 


SaKe I I(t) = Bowe + UES 
eG el P eek Hoh hans peg eat /k0 
dt 
t — _ 
(pe ae ar ipa 
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These curves can then be analyzed by simply looking at the 

regions where simple exponential behaviour is evident or by 
using the multi-exponential fitting method. Although these 
techniques appear promising at first sight, in practice they 
do not give more reliable values than the multi-exponential 


curve fitting method discussed previously. This is because 
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of the problems associated with the numerical differen- 
tration of statistical?ddatas 

Another possible technique is to measure the decay of 
alignment of an excited state (BCS 72). This method assumes 
that transfer of alignment by cascades is small. A Similar 
method, which is applicable even if cascade alignment does 
occur, uses the decay of g-value quantum beats ({LC 72). 

Both these methods use only 1-10% of the available intensity 
and are of limited applicability. Finally, a coincidence 
counting experiment (MS 73) has also been tried, but the 
counting rate was found to be much too low for accurate 
lifetime determination. 

A potentially useful technique was suggested by 
iedsoCurtis,andshisxcolleagues an) 1970« (CBBI70,,€CBE, 71). In 
this technique the decay curves for the primary transition 
and all its direct cascades are recorded and then analyzed 
Simultaneously. In principle, this method can completely 
account for the effects of cascading on the value which is 
determined for the lifetime of the primary level. This 


method of analysis is the topic of the rest of this work. 
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CHAPTER Bat 
THE ANDC METHOD 

The application of the ANDC! theory to be derived in 
section 2.1 will, in principle, solve the problem of cascade 
effects in radiative lifetime measurements. However, the 
presence of experimental errors (i.e., statistical 
fluctuations of photon counts) and the finite amount of data 
means that there will be some error involved in the 
determination of the lifetime. Thus the problem is not only 
how to apply ANDC analysis to best advantage, but also to 
provide reliable error limits on the lifetimes obtained. 
Section 2.2 describes the calculation of the basic ANDC 
variables and associated errors, and also provides an 
introductory discussion of the relative merits of the 


different methods of caiculating these variables. 


Consider the generalized decay scheme in Fig. 2.1, 
where a primary level (p) decays into a state (m) and the 
primary level has a number of higher levels decaying into 
it. These higher levels are also subject to cascading. By 
including ail direct transitions (n+p) into the primary 
level, the population rate of change of the primary level is 
found to be 


1 Acronym for Arbitrarily Normalized Decay Curves. The 
reason for this name will be made apparent in section 2.1. 
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aN, {t) 
2 IN, (HAL, a ey pen - (2.1) 


Fig. 2.1 Generalized Decay Scheme. 


Noting that the observed intensity for the transition i+j is 
rig i (Aas where Ei is the relative detection 
efficiency! of the experimental apparatus for the 


transition, the following expressions can be substituted 


into equation 2.1 3 


eee T(t) an (t) . a©. /at 


BE. A ELA ela EA 
pm pm np np pm pm 


This gives the result 


Meal Ob : if em I 
EL eA 
pm pm 


which can be written as: 


1 The detection efficiency is determined by the optical 
properties of the experimental apparatus used and depends on 
the wavelength of the transition. 
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dI EA 
Tp —= = }Enl_ - 1,,> where Ep = SR PE . (2,2) 
dt lia ve 
np 


Integration of equation 2.2 with respect to time between 


initial time, te and final time, te yields: 
BD } = f ; 
ae ont e) tom ca = ~ {Todt + bend tpt ; 


Ox ray i ere here (253) 


where P = | ber gaor sk area under primary curve, 


Cc =i Gee vevat = area under cascade curve 
and Lae Onin - ice = change in primary 
intensity 


This last equation is linear in the variables P/, and 
C/A, whose values depend on the choice of t. and t. (see 
Fig. 2.2), with constant coefficients 1 and €. respectivelyt 
The constant term is the primary lifetime, which can now 
be obtained by performing a linear fit using the variables 
Pf and Cy. If no errors were present in the data, a ,_ 
simple linear fitting procedure could be used. However, as 
mentioned previously, the presence of experimental errors 
requires the use of special fitting procedures to account 


eee ee es ee EY 


1 Note that since the €, contain the relative detection 
efficiencies of the transitions, the absolute normalization 
of the decay curves is irrelevent in determining the 
lifetime. 
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properly for these errors. 


It is possible to write equation 2.3 in the form 


N= 


alt 


= : = clare (2064) 


This equation is not as simple as 2.3, since it has an 
additional variable term. However, it does have certain 
advantages due to error considerations, as is shown in the 


next section. 


222 Practical Application of ANDC 

The basic quantities (P,C,4A), which are the areas under 
the decay curves and the change in the primary intensity, 
are most easiiy calculated directly from the signal data. 
To do this, the cascade and primary intensities must be 
measured at the same distances from the foil. If this is 
not the case, either interpolation between data points or a 
smooth fit to the cascade decay curve can be used to 
calculate the cascade intensity at distances determined by 
the primary data.! Linear interpolation between distance 
coordinates is preferable to fitting the curve for the 
reasons given below and in Chapter III. 

The principal objection to using a fitted curve, 
1 The ANDC eae the case of one direct cascade 
presented later in this work allows the six parameters of a 


HOMER multi-exponential fit to be used to "reconstruct" the 
cascade data set. 
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particularly for the primary, ie that any subsequent error 
calculation in the ANDC analysis may be of doubtful 
validity. The reason for this is that the estimated errors 
in the fit parameters usually reflect how good the fit is as 
a whole, even though the fit may be poor (or excellent) ina 
particular portion of the decay curve. In other words, 
Since the true decay curve is a complicated sum of 
exponential terms, which is difficult to fit even without 
statistical fluctuations, any condensation of the 
information into a few parameters will generally lead toa 
loss of information about particular parts of the curve. 
This loss can be particularly important when calculating the 
intensity difference, A. 


The areas are calculated using the trapezoidal rule: 


f-1 
Cit aired) acts apsiay ce) | 
P= ) boast nah cot ea mawienre Y= intensity. 
— 2 
Je 


Thus the variance! in the area is given by 


f-1 2 
Yeeros Ein ares ee 
ea) rd fa ae 
— 4 
je 


Since there is no correlation between the intensities at 
different times and a ea due to Poisson statistics. 


The difference, A ¥, $Y has a very large relative 


aha s wate 


error when compared to that of Pj rf or Cyr, since V (4) =, +¥_-- 


by ae 


1 Notation: variance V(X)=<(dX)2>. 
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For exampie, if Y,=9000 and Y-=8000, 


L 
S(P) /2 (9000+8000) *7/17,000 = .008 , whereas 


FE 
S{A)7A (9000+8000) */(9000-8000) Sebo Uris 


This is especially true when A is small. To minimize this 
error, the intervals toSte(referred:tto asipanel) sizes) are 
made as large as possible and a single exponential fit is 
performed upon the five point section of the decay curve 
centered on the desired point in order to obtain Y,; or Yee 

This fit is necessary in the tail of the decay curve, 
Since it can happen that the measured A is less than zero 
due to the statistical fluctuations. The fit also reduces 
the error in A since five points are being used to obtain 
the intensity at one point. Thus we can divide the variance 
by the number of degrees of freedom (5~2=3). 

There are two possible ways to proceed from this point. 
Considering the simplest case of one direct cascade and 
referring to equations 2.3 and 2.4, one can either fit 
ppanES. Of !the form APC) OR (P/N, C/A). “The Latter’ choice 
looks simpler in that it requires only a two dimensional 
linear fit, whereas the other requires a three-dimensional 
fit. However, dividing both areas by A not only introduces 


large error to both variables, but also causes almost 100% 
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correlation between the errors of the two variables.! 

Assuming a linear relationship between X and Y, it is 
possible to write an expression for the probability of 
obtaining the set of data points (X,Y) and associated errors 
(<dX>,<dY>,<dxdY>). This probability will be a function of 
two parameters, the slope and the intercept of the line. 
Maximization cf this probability with respect to the two 
parameters will yield equations for them which can be solved 
numerically. The errors in the two parameters can be 
estimated by a similar procedure (CRRW 72). This procedure 
results ina satisfactory lifetime and associated standard 
deviation in caseS where it is applicable, that is, when 
there is only one significant direct cascade present. 
Examples of its application are given in subsequent 
chapters, where this is referred to as the (2-D) method of 
analysis. 

In the case of only one direct cascade, one can 
eliminate much of the correlation between the variables of 
the two dimensional fit by writing the ANDC equation in the 


forn 


1 There is also a small correlation between P and A, but it 
has a negligible effect on the results. 
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This reduces the correlation of errors of X and Y and 
also reduces the error in X, since the largest source of 
error, A, is now only used to calculate Y and not X. 
However, tests on real data showed no significant difference 
the lifetimes and errors obtained using this equation! 
instead of the (2-D) method, when the method is applicable. 
This occured despite the fact that the correlation 
coefficient, R=< (dX) (dY¥)>/(S(X)S(Y)), for each panel was 
reduced from about 0.95 to 0.1 in going from the (2-D) to 
the (2-A) method. This is a good indication that the error 
calculation and weighting procedure of the linear fitting 
program is valid, even when the correlation between the 
variables is very large. A test was done on a case which 
had two direct cascades, but the analysis was made including 
only one. Both methods of two dimensional analysis produced 
unacceptable (put different) results. The alternate two 
dimensional method may be of use when the error estimation 
procedure of fitting routine has trouble converging. 

In the case of the three variables (A,P,C), there is no 
significant correlation among the variables. Since the 


errors in the areas are very small (<<10%) compared to the 


1 Referred to as the (2-A) method later on in this work. 
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error in A, the multi-linear fit can be weighted using only 
the variance of A.4 The fitting procedure is more 
Straightforward than the two-dimensional case above, but the 
errors cannot be accounted for rigorously. The obvious 
advantage to this procedure is that it caa easily be 
extended to cases where more than one direct cascade is 
present. 

There are also two distinct ways of partitioning the 
decay curve into several panels. The panels can be 
overlapping or non-overlapping. The latter choice severly 
crestricts the size and number of panels (i.e., data points), 
but produces data points which are completely independent of 
each other. Overlapping panels will resuit in some 
correlation between data points, which will make the fitting 
procedure and error estimates obtained from it less 
rigorous. 

The next two chapters will be devoted to testing these 
possible variations of the ANDC method on synthetic data 


(Chapter III) and actual experimental data (Chapter IV). 


i This is very convenient, as there’ is no statistically 
valid procedure for performing a multi-linear fit if the 
errors in the variables are comparable. The computer 
program used is REGRES ,which is taken from Bevington (Be 
69). 
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CHAPTER III 
ANDC ANALYISIS OF SYNTHETIC DATA 

This chapter establishes the usefulness of the ANDC 
method of analysis in situations where significant cascading 
is presente Conventional methods of analysis may give poor 
Or misleading results in these cases and ANDC can serve as a 
valuable check on the results. 

In order to assess the accuracy and reliability of the 
different ways of applying the ANDC method described in the 
previous chapter, a few appropriate situations where the 
lifetimes are exactly known are necessary. This is most 
easily accomplished using artificially constructed decay 
curves. Comparisons of the ANDC results with the known 
values and the results obtained using conventional curve 
fitting techniques reveal that ANDC is an important tool in 


the analysis of cascaded situations. 


The model decay schemes used are shown in Figs. 3.1- 
3.3. The initiai populations of the levels were chosen to 
be proportional to (21+1)/n3, which previous studies 
indicate is approximately the case in many physical 
Situations.!1 Branching ratios and lifetimes were chosen to 
represent some typical situations where cascading effects 
are important. The analytic expressions for the intensities 
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Fig. 3.1 Synthetic Data: case 1, 
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Hig. 3.2.oynthnetic data: case. 2. 
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Fig. 3.3 Synthetic data: case 3. 
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were obtained from N,;(0), Ai j and T; using the general 
formula given in Appendix I and assuming Ty 5=AigNi (t)- 

Each data set has sixty points and all sets 
corresponding to the same case have the same set of distance 
coordinates. The program used in the analysis assumes that 
the data are subject to Poisson statistics, and it 
calculates errors and weights points accordingly. This 
fluctuation is present in all real data, so an atitempt is 
made to make the synthetic data fluctuate about the "true" 
value according to Poisson statistics by using the FUZZ 


subroutine of the program HOMER (IL 74). 


3.2 Analysis 

The analysis was performed in several different ways in 
order to determine which would give the most accurate 
results. For the two cases (Figs. A aba: where only one 
direct cascade was present, both the two-dimensional (2-D) 
and multi-dimensional (M-D) methods were used. The (2-D) 
method was not used for the case shown in Fig. 3.2, since 
the presence of a second direct cascade would render the 
results inconclusive. 

The data points were partitioned into panels three 
different ways, with a varying degree of overlap between the 


panels. The number of ANDC points! was kept at about 


1 The term "ANDC point" refers to the calculated quantities 
(P,C,A) associated with a particular panel(cf. Fig 2.2). 
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fourteen for each partitioning. Table 3.1 gives the three 
different partitionings used, (A,B,C), where the numbers 
define the endpoints of the panels and refer to the data 


points in order of increasing distance from the foil. 


Table 3.1 Panel divisions for synthetic data analysis. 
A paneling contains 15 ANDC points: 
fay OeOPI5.20825-30 85040 "45050859532 37420 
Ome 20.2) S083 5.550 045 42: 752) 56> 608 Bb. 40) 60 
B paneling contains 14 ANDC points: 
toes COSISh 479 2482984538 Fond 1920030050 
feo 2 16 25 26) :33 37745590 Su 30-4060 


C paneling contains 14 ANDC points: 


Pp 
nme OVE OLAS Tie 2itez Se29"sses7041746) 51 56 
Demoneve., 16) 207.24. 28 182360740) 45) 501555: 60 


The first partitioning, A, has a large degree of 
overlap between successive panels, which allows for large 
panels each having 10-15 points. The second partitioning, 
B, has non-overlapping panels for all but the last three 
ANDC points, which overlap most of the preceeding panels. 
The last partitioning, C, provides no overlap at all, but 
conseguently gives smaller panel sizes, typically having 
about 5 points per panel. 

The program yields the primary lifetime in each case, 
plus quantities which would be the branching ratio in the 
(2-D) analysis or the transition probability in the (M-D) 


analysis, provided that the detection efficiencies for the 
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primary and cascade were the same (cf. eqns, 2.3,2.4) The 
estimated errors in these quantities are also obtained, 
along with the value of the reduced chi-squared for the 
fits! 

An important question now arises. How is it possible 
to decide when a particular analysis is unacceptable? An 
unacceptable fit may occur due to poor data or when the 
analysis performed 1s inconsistent with the actual physical 
Situation, e.g., the analysis is performed assuming only one 
direct cascade when two are really present. 

One indication of the validity of the analysis is the 
value of the reduced chi-squared (y?). The probability that 
the fit obtained is valid is small when the chi-squared is 
Significantly greater than one. 

In the multi-dimensional analysis, the fit allows for 


an “extra"™ constant term, K, 


which should be equal to zero. If the value obtained for K 
differs significantly from zero (considering the estimated 
error in K and the typical value of A used in the analysis), 
then the analysis should be rejected. Also, the multi- 
dimensional analysis assumes that the error in the area is 


1 See Appendix IIi for a more detailed description of the 
program output. 
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negligible compared with that in A. The validity of this 
assumption should also be checked in each analysis. In the 
two-dimensional analysis, the plot of P/A versus C/A should 
consist of points randomly scattered about a straight line. 
The presence of a blend or additional direct cascade will 
generally result in a curved rather than a straight line 
plot, thus giving another criterion for rejection of an 
analysis. The sign of the transition probability (modified 
by efficiency factors) should obviously be positive, giving 


a further check on the validity of the fit. 


The resuits of the analyses are shown in Tables 3.2-3.4 
and reveal several important features. Most of the values 
obtained for the lifetimes, transition probabilities and 
branching ratios fall well within the estimated single 
standard deviations and all values fall within two standard 
deviations of the true value. It is apparent from Tables 
3.2-3.4 that the standard deviation is dependent on the 
average panel size; i.e., the larger the typical panel size, 
the smaller the standard deviation. The results also 
indicate that the estimated standard deviations on the 
lifetime values are consistent with both the deviations of 
the individual values from the true result and the scatter 
among the four values obtained for each individual 
partitioning pattern (A,B,or C) used. The degree of overlap 


given by the A partitioning appears to be close to the 


pay 
ui? a 
Lit 3 
§ i 
e 
os 
? 

I i 
ii ? 45 
4 . «- 


26 Castn 


a 


= @ » oe P 
sing b iS Ee 


> 
i 


2 hetar * 


euarey. B 


1 inte s Poaee 2 SORE or ends 


»y SON > 


Ke 


AT oot pig eb, Storrs te say sali 


AES 262 be ute: » sou lev gue? sea paows _ 


 oheie vader askew beaerann Che ig 
4 Gi fodder oa) vets — 


Jolag. oko ‘seine Lene 


a Date save "" nay fates, on or toMy sa 
j p SESSA a he ee. * ae aco ae 
ne Hg tie Oe sesnaes yites 


pred” MQoee, 4696 re’ coven 
o (eceefen ot? 46 tpi ai 


R2AS GE teE RS ata 
enebeire + a hee 

the shkdanga’ Dea vee Lagve 
zoe onkre ret sikh yak hou: a ; di 

odd argnew RM bee) wes pa 
beth oe ibl spot Pee 6 aedeheeb ‘Bs | 
lee Han tak tn eno 
Loic ivedi tub anee ee Sets BoE 


aa Pits OVS wettel DAS LE rnsiaeee 


f Logos hy@aurias oat. tana ® Sapkhad: 


nat es 
ind otis seedeicns 34s conte’ see 
> t sW7 mot? boyde® Leber eee a 


beep (2 TOG eh) Saas, vitontadiaag 
(F 2: boqgs  BRRAGESTS ae 4 ‘ oe pa vow Lp 


Table 3.2 Synthetic data case 1: ANDC lifetimes. 


PERM ii 
AA 6.90 
AB Ge T1 
BA ha t6 
BB 7. 14 
AA 7.61 
AB Fe 65 
BA res Ws 
BB Te 29 
AA Ya ¥2 
AB 7243 
BA al 
BB hat2 
AA 6.68 
AB 6.64 
BA 6.78 
BB 62:65 
AA Tw 2g 
AB 4229 
BA 6.76 
BB 6.19 
AA 6.49 
AB 6.243 
BA 6.66 
BB 62172 


SYN#1 (7T=7. 14, A=0 14, BR=. 99) 


O 
0.36 
0.36 
0. 38 
0.39 


0.2 68 
0.69 
0.67 
0.70 


1.01 
1.05 
1.02 
1. 16 


Oe o2 
0.53 
0.50 
Ge. 32 


0.73 
0.75 
0.72 
0.74 


0.94 
0.94 
1.09 
ASag 


A 
0.15 
0. 16 
0.14 
0.14 


0. 13 
0.13 
0.14 
0. 14 


0. 14 
0.13 
0.14 
0.13 


BR 
1.04 
1.06 
1.03 
1.05 


0.98 
0.98 
0. 92 
0.94 


71.04 
1.06 
12 03 
1x03 


6A 
0.01 
0.01 
0.01 
0.01 


0.02 
0.02 
0.02 
0.02 


0.03 
0.03 
0.03 
0.03 


6 BR 
0.05 
0.05 
0.05 
0.05 


0.07 
0.07 
0.07 
0.08 


0.09 
0.09 
0.10 
Oe 11 


x2 
0.97 
1.00 
1s37 
T39 


0.71 
0.88 
0.87 
1222 


0.53 
0.62 
0.74 
1.07 


0.76 
0.82 
1.01 
eZ 9 


0.79 
0.83 
0.63 
0.84 


0.54 
0.59 
0.44 
0.63 


(M-D) 
A PANELS 
extensive 
overlap 


("-D) 
B PANELS 
MODERATE 
overlap 


("-D) 
C PANELS 

no 
overlap 


(2-D) 
A PANELS 


(2-D) 
-B PANELS 


(2-D) 
C PANELS 


optimum, giving the lowest estimated standard deviations 


with reasonable values for the reduced chi-squared; 


ise., the panels are sufficiently large to give small 


errors, while maintaining a reasonable degree of 


independence between adjacent panels. 


Another important feature 1s the pattern followed by 
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Table 3.3 Synthetic data case 2: ANDC lifetimes. 


SYN#2 (tT=7.14,A=. 14, BR=.99) 


PERM T Oo A oA x2 
AAA 7203 0.53 0275 0802 Tat «| 
NAB 8660.87 0.52 O475 0007 AGE ii 
ABA 6.76 0.53 ONTS 0202 1% | 
BEB ~6s.70 0.52 On) 1S .de02 ieesaie (M-D) 
BAA 7.280 0.61 O212 0802 -63 {| A PANELS 
BAB 7.265 0.60 OF16 0802 s52 4 
BBB 7245 0.60 0213 0202 25h | 
BBA 7.50 0.60 04138 0402 258 «| 
AAA 9.16 3.30 0.08 0.06 Pe aay 
: a we $ | (M=D) 
- : AO ws : | B PANELS 
J 5 are k 
REJECTED DUE TO LARGE ERRORS | 
AND LARGE CONSTANT TERM 
ABA’ | Ba 25 O96 0.14 0.03 Bic ee 
RAB 7.16 0.96 Ofte 0.203 -68 | 
ABA 6.78 0.95 0.15 0.03 WHO | 
ABB 6275 0.96 O25 0203 Syl at (M-D) 
BAA 8259 184i GATO 0203 =36 § C PANELS 
BAB 8.49 1.40 0.10 0.03 538 f 
BBB 8.20 1.50 ORT? 0503 oth | 
BBA 8223 1248 OF oO fot =o 


the permutations of primaries and cascades. The value of 
the lifetime obtained depends, in most cases, mainly on the 
primary and only slightly on the cascade matched with it. 
This is not unexpected since the primary contributes two 
variables (P and A) whereas a particular cascade affects the 
result through only one variable (C). Furthermore, the 
major source of uncertainty lies in , which is derived fron 
the primary decay only. 


In order to obtain the best estimate of the lifetime, 
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Table 3.4 Synthetic data case 3: ANDC lifetimes. 


PERM 45 
AA 5-40 
AB 5 43 
BA 5.44 
BB 52 46 
AA Sa bs, 
A3 Dent 
BA 5. 46 
BB Spee 
AA De (3 
AB 52 14 
BA 5.42 
BB 5242 
AA 5339 
AB 5.43 
BA 5. 43 
BB 5245 
AA Dat 
AB Dersz 
BA 5-50 
BB ee a 
AA 5.43 
AB 5. 46 
BA Pye) 
BB 25 57) 


SYN#3 (T=5.40,A=.125,BR=. 675) 


Q! 
O. 11 
O. 11 
On 12 
0.212 


0.23 
0.23 
0.23 
0. 23 


0.31 
Us 1 
05.33 
0.33 


Ont 
0.11 
0.11 
Osi 


DNC 
DNC 
0.17 
0.17 


O29 
0. 30 
a Ps 
DNC 


OA OA 
2 124 50.06 
2-124 .006 
2122 .007 
aA2L ed 


Pe Je > Nestea 8) 
2139 2012 
eh et 
eS 120 0A 


2M. 08 
e120 2018 
2119 .018 
leno oreo 


BR BR 
36594. 022 
-660 .024 
2658 .024 
3063/6025 


mOI0YDNE 


-697 DNC 
639 .040 
-645 .040 


BO Oe ie 1 Ow 
500 4.106 
e594 2110 
2994 DNC 


~19 


(M-D) 
A PANELS 


(M—D) 
B PANELS 


(M-D) 
C PANELS 


(2-D) 
A PANELS 


(2~D) 
B PANELS 


(2-D) 
C PANELS 


DNC--> error analysis routine did not converge. 


some method of averaging the results of the different 


analyses performed is required. 


Of the results obtained, 
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only the B paneling on SYN#2 (Tabie 3.2) had to be rejected 


according to the criteria described above. 
were rejected due to the very large estimated standard 


deviations (~30%) and also the occurence of negative 
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transition probabilities. For the remaining analyses, the 
X* of each result was included by multiplying the estimated 
standard deviation in the lifetime by X2 when X2>1. Using 
these modified error estimates, the final weighted average, 


ct, and its estimated standard deviation, o , were calculated 


using 
ite Til 
Io gt ab 


The results for all three case studies were reduced 
uSing this method and the final values obtained are given in 


Table 3.5. The ({2-D) and (M-D) methods were averaged 


fable 3.5 Lifetimes of synthetic data: summary. 


HOMER 2-D 

CESS 2-D TROY M+D M-D 

#1 6. 80 Se 7217 6.98 

T=72.14 £522 +2 to20 +.29 

3% 25% 3% 4% 

#2 Nei Ae Bett Pi Nas 
t=7.244 +. 39 +245 
5% 7% 

#3 5.44 Sy 5.40 52 42 

T=5. 40 +.04 te 16 £5°10 +.08 

1% 3% 2% 2% 


separately and together where applicable. The results of 


the conventional method of analysis are also included for 
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comparison. The programs HOMER and TROY! were used for the 
multi-exponential analysis. TROY is a program which fits up 
to six exponentials plus a constant term to a decay curve. 
Unlike HOMER, it allows any number of the thirteen possible 
parameters to be fixed. This can be useful for finding the 
primary lifetime when the cascade lifetimes are known. 

Using the cascade and the primary lifetimes as determined by 
HOMER and/or ANDC as initial estimates for the parameters of 
the multi-exponential fit, it 1S sometimes possible to 
confirm the ANDC lifetime estimate. This is done by holding 
any combination of these lifetimes fixed to obtain a fit 
consistent with all the "known" lifetimes. 

Table 3.5 shows that the ANDC analysis provides quite 
Satisfactory resuits, as regards both accuracy and 
ELrecasion,@invall three’cases. It 1s worth noting that, as 
one might expect, the percentage error in the lifetimes 
increases as the decay scheme becomes more involved. 

The method of ANDC analysis described in this chapter 
requires that both the primary and the cascade data have the 
same set of distance coordinates at which the intensity is 
measured. Unfortunately, not all experimental data is in 
this form, so the original cascade data set may have to be 
modified. This can be done by either simple linear 
interpolation between points or by use of a multi- 
exponcntialerit to the cascade. “it. is possible for the 
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latter choice to significantly reduce the reliability of the 
ANDC results for the reasons explained in Chapter II. This 
is demonstrated by the result of doing the previous ANDC 

analysis with fitted curves rather than the original cascade 


data. The result is shown in Table 3.6. The lifetime 


Table 3.6 ANDC lifetimes obtained with original and fitted 
cascade data. 


CASE ORIGINAL FILET ED 
Gis 90) G36 6094 37 | 
# 1 Go: 16,36 6250 «37 ( A 
7.14 hetl6 tat3 3 tO) cele] {| PANELS 
72.14 239 Cisil 2 sia | 
Cal2 le Gai tz «190 { Ss 
7243 1.1 G0 Jee 0 { PANELS 
iD atiths. Ae O 6.59 1.0 1 
TMZ WeNierZ 6) 516)7 te 97 i 
AV ERAGE: 7.20140.2 6.784+0.1 
(3%) (1.5%) 
5-40 11 5243 2.11 { 

# 3 5243 211 547 211 | A 
5. 40 5044 212 Sak] bei hZ } PANELS 
5at6. 2 12 5049 212 1 
epicure od D8 INS ee | 

5214 2.31 SU ve Su | C 
5042 . 33 5e44 . 33 | PANELS 
5.242 «33 59933 | 
AV ERAGE: 5.4124. 07 5244+. 08 
(1.4%) (1.6%) 


obtained for Case 1 from the analysis using the fitted 


cascade data set is misleading in that the error estimate is 
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much too small. However, the analogous resit for Case 3 is 
essentially identical to the result obtained using the 
original cascade decay data. 

Linear interpolation is more reliable since it cannot 
produce any systematic deviation from the actual 
experimental data (provided that the data points are taken 
at reasonable intervals). The linear interpolation will 
only affect the exact manifestation of the Poisson 
statistics of the data by reducing the random fluctuations, 
ieee, the cascade decay curve will be partially smoothed in 


the process of interpolation. 


3-4 Summary 

On the basis of the preceding analyses of synthetic 
data, several rules can now be given for the application of 
the ANDC method to real data. Firstiy, it is more important 
to maximize the guality of the data (counting statistics, 
Number of runs etc.) for the primary decay rather than for 
the cascade decays, since the calculated lifetime depends 
principally on the primary data set. Also, where posssible, 
the data should be taken so that the distance coordinates 
are the same for the primary and the cascade(s). If this 
last condition cannot be fulfilled, then it is preferable to 
use linear interpolation rather than a multi-exponential fit 
for the cascade. 

During the actual analysis, the paneling should be done 


in at least two different ways, one with no overlap and one 
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with moderate (50%) overlap between adjacent panels, 
€ege, Table 3.1, C and A panels respectively. The results 
of both panelings should be consistent, even though the 
overlap analysis will probably give the better results. It 
is also useful to check the results by using both the (2-D) 
and (M-D) analyses (where possible). 

The following quantities can indicate that a particular 


analysis should be cre jected. 


i) reduced xX@ --reject analysis if too iarge, 


1i) curvature of graph (2-D analysis only) --reject if 
Signiiicant curvature present, 


foi) validity" (M—D only) ==reject if these quantities. <1 
for too many points in the analysis, 
ev) constant K’ (M-D only) -- reject if too large and/or 
estimated error in K is too small (considering the 
typical value of A); e.g., if Av200 and K=8024E; one would 
definitely reject this if E<10 and would possibly accept 
it (depending upon the other criteria above) if #240. 


v) sign of transition probability -- reject if negative. 


It 1s important to bear in mind that these rejection 
Criteria are not infallible. Ityis aways possible for 
circumstances to combine to give a physical situation which 
yields an “acceptabie" result from an analysis which does 
not correspond to the physical situation in the experiment. 
For example, a significant direct cascade may have been left 
out of the analysis, yet the analysis may still be 


"acceptable" according to the criteria given above, thus 
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yielding a misleading value for the lifetime. An important 
safeguard in this regard is a study of the energy level 
scheme and the cascade scheme relevant to each particular 
primary iifetime measurement to insure that no significant 
cascade transition has been overlooked. No amount of 
mathematical manoeuvering on the computer can compensate for 
an oversight made in the laboratory. 

In Chapter IV, examples will be presented showing the 
application of these general validity criteria to ANDC 


analyses performed on actual beam-foil decay curve data. 
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CHAPTER IV 
ANDC ANALYSIS OF REAL DATA 

This chapter provides examples of the application of 
ANDC analysis as described in the previous chapter to some 
real dates. The examples include cases which were intended 
to be used in conjunction with ANDC analysis in that the 
primary and cascade distance coordinates are matched, and 
also cases where linear interpolation and/or multi- 
exponential fits had to be used for the cascade(s). 
However, before these examples can be discussed, one must 
first consider the possible effects of the instrumental time 


averaging that occurs in the collection of actual data. 


4.1 Dime Averaging Effects 

Experimentally, one does not measure the intensity at a 
particular distance, d,, from the foil; rather one measures 
an average intensity over a segment of the beam centered on 
the point, d,- Fig. 4.1 is a schematic of the experimental 
Situation when taking measurements below 2000A and it shows 
that the averaging time, 6, depends on the mask size, m, 
slit isaZe, s, and the beam velocity, v. Mask sizes in this 
laboratory usually range from 5-52mm, slit sizes from 100- 
800m, with typical velocities in the 1-2 mm/ns range and 
d,/d, 0-06. Note that, due to this averaging time, the 
first meaningful data point after the foil occurs at time 
t=6/2 (unless one is able to include the experimental window 


function in the data analysis, THB 77)< .ethis results in 
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Fig. 4.1 Optical configuration for wavelengths <20002,. 
(not to scale) 
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values for 6 ranging from about 0.2 to 4 nse. Fig. 4.2 shows 
the analogous situation when the wavelength is greater than 
2000A. In this case 6=S/v and s<800 m, so this case is 
included in the previous one. 

Frequently, 6 is made larger than is desirable, sacri- 
ficing good time resolution in order to obtain a suffic- 
iently high counting rate. Thus when measuring lifetimes in 
the range less than 1 ns, 6 is often as large as the value 
of the lifetime being measured. The question now arises as 
to whether this will affect the value obtained for the 
lifetime and, if so, what value must § have relative to the 
lifetime so that the effect is significant? 

As an example, take a two component decay curve with 


"ideal" intensity 


Integration over the segment, §, yields the measured 


intensity 
@ shane 26205072. Ot ge 
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+B sh (a26/2) e725 /2 aeons 
Q2 2 
where the factors of nono f2 and a7t2/2 are included so 


that the first valid data point corresponds to t=0. These 
two equations show that the shape of the decay curve will 
change depending on the values a16 and a25. Taking S, the 


ratio of the first and second coefficients, as a measure of 
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the shape of the curve, 
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The shape of the "ideal" curve is S;j=A/B. 
A quantity which will reflect how the measured shape 


deviates from the ideal case is 


The value of R for some values of a;5 anda2dé are given 


in Table 4.1. It iS apparent that in most cases the shape 


Table 4.1 Values of the Shape Factor, R, for some a.6« 
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multi~exponential fitting in order to obtain the lifetime 
Since the various exponential terms are unchanged. But what 
happens to the ANDC analysis? 

ANDC analysis states that, in the case of one direct 


cascade, the lifetime of the primary level is given by 


If A is such that one or both decay curves change 
Shape, then the quantities (P,C,and A) become (P',C*, and 
NE) Pisiiitct hen +t cue’ tthat 
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If not, how much deviation from the true value of the 
lifetime can be expected? 

An analytic calculation was done which suprisingly 
enough shows that the value of 6 has absolutely no effect on 
either the value of T or € obtained from ANDC analysis. The 
details of this computation can be found in Appendix Iil. 

Of course, in practice one cannot make 6 arbitrarily 
large for basically the same reason as one cannot do so when 
using a multi-exponential fit, the reason being that the 
information from the component with lifetimes much less than 
§ would be lost in the larger random fluctuations of the 
data, resulting in very large uncertainties in the values 


computed for then. 
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TO summarize this section, one can conclude that there 
is no special difficulty in applying ANDC analysis to decay 
curves which have been 6 SOU affected by time 
averaging. Of course, just as for muiti-exponential curve 
fitting, the uncertainty in the primary lifetime will 
increase as the averaging time, 6, is increased. 
Furthermore, for the ANDC analysis to remain valid, it is 
necessary that the same averaging time be used for both the 


primary and the cascade(s) in a given ANDC data set. 


4.2 Fluorine Ii!: 3p3D Level 


The relevant cascade decay scheme involving the 3p2D 


level is shown in Fig. 4.3. The ANDC analysis was performed 


Fig. 4.3 Decay scheme of F II - °D level (numbers in {} 
are approximate intensities). 
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1 All the data used for the analyses in this chapter were 
obtained at the Radiation Research Laboratory of the 
University of Alberta. 
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using the transition (3973A) to the 3s3P level as primary 
and the line at 3588A as the only direct cascade. The line 
at 3125A was not detectable, which indicates that cascading 
effects due to this transition were small. The direct 
cascade at 3544A was severly blended with the transition 
3d3D-5p3P and 3s3D-3p3P and therefore could not be included 
in the ANDC analysis. It was, however, weak relative to the 
3588A line and therefore should not be a major contributor 
to repopulation of the 3p3D level. 

Three runs were done on both the primary and the 
cascade, giving nine permutations of primary with cascade. 
The beam veiocity after the foil was 2.08mm/ns (+2%) with 
the beam energy at about 1.0 MeV. There were between 40 to 
48 data points in each of the primary data sets with peak 
counts of around 4400. The background count rate was 
Similar for all three primaries (~4 per sec.), but the 
average coilection time per point was a factor of two larger 
for primary #1 (P#1). The cascade data sets have 38 data 
points each and peak counts between 3000 and 4000. The last 
data point taken, for both primary and cascade, corresponds 
to a time since excitation of about 40 ns. 

Two types of partitioning were used for the ANDC 
analysis, as indicated in Table 4.2, which also contains the 
results of the analysis. These data were not taken with 
ANDC in mind, so the cascade decays were replaced with 


fitted curves in the analysis. Both multi-exponential 
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Table 4.2 Results of ANDC Analysis on F II Data. 
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order to compare the results of each when used on actual 
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little difference between the two methods of fitting the 
cascade, the linear interpolation fits providing only 
Slightly longer lifetime values and slightly lower reduced 
chi-squared values than the multi-exponential HOMER fits. 

The two dimensional analysis(2-D) was the principal 
method used, but for comparison part of the analysis was 
done using the multi-dimensional analysis(M-D). The (M-D) 
analysis produces consistently larger values for the 
lifetimes (Fig. 4.2) than the (2-D) analysis, although the 
two agree to within the error estimates. 

All the permutations of cascades with primary #1 were 
rejected due to a large reduced chi-sguared. This is the 
result of the poor quality data of P#1, as attested to by 
the longer accumulation times. All the other permutations 
were acceptable according to the rejection criteria 
described in the previous chapter. Figs. 4.4 and 4.5 are a 
typical example of the results produced by the computer 
analysis! for this case. 

The overall average lifetime was first calculated using 
only the starred(*) values in Table 4.2, i.se., using only 
the best fit for each permutation. This gave a value for 
the“lifetimne of 6.9640.2 ns. The ré@sults of all the valid 
analyses were then included in the average and produced 
essentially the same value 6.994#0.2 ns. 


This latter value of the lifetime was then used, along 


1 See Appendix II for a description of the computer output. 
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With a cascade lifetime of 4.49 ns obtained from a HOMER 


Single exponential fit, as initial estimates for a TROY fit, 


holding the cascade lifetime fixed. 


Good fits were obtained 


for primaries #2 and #3 with reduced chi-squared values less 


than 0. 


88.6 


This gave a value for the primary lifetime of 


72024+0.22 ns, in excellent agreement with the ANDC result. 


Fige 4.4 Sample analysis(2-D) for 
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This 1s an example of the use of two direct cascades in 
ANDC analysis. The relevant part of the decay scheme is 
Shown in Fig. 4.6, and Table 4.3 lists some of the 


characteristics of the data obtained. 


Fig. 4.6 Decay scheme for I VII - *Pay. level (numbers 
in {} are approximate intensities) . 
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HOMER fits had to be used for the cascade since all the 


Table uss Characteristics of 1 Vil /data-. 


PEAK ACCUMULATION INTEGRATION 
(A) #RUNS COUNTS TFIMECPERS POINT) LENGTH(NS) 
PRIMARY 834 2 3700 18 SEC Onis 3 
CASCADE 640 Ps 2000 TSEC 027s 
CASCADE 464 1 2050 355556 Os66 
Beam energy 1.3 MeV Beam velocity=1.537 mm/ns +1% 


Last data point used in analysis occurs at t=4.4 ns. 
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curves were measured at different distance coordinates. The 
E03: (9-- P3795 transition at 649A (Fig. 4.6) was not included 
in the analysis since the line was severly blended with 
another stronger line and is also expected to be weak 
relative to the EG ly boeing) cascade. Its decay curve should 
also have a very Similar shape to that of the Ds hoe Be 
Pransitione 

The transition at 464A (5p-6s) was included in the 
analysis even though an examination of the spectrum 
indicates that it probably has a significant blend with the 
4€-5g transition. An attempt was made to do the analysis 
Without this line, but all the results had unacceptably 
large reduced chi-squared values, even though the the 
lifetimes obtained were close to the those obtained when the 
464A line was included in the analysis (this is not 
unreasonable since this cascade had half the intensity of 
the 5p-5d cascade). The lifetime of the blend may be such 
that thershape of’ this’ cascade* 1s" not Significantly changed. 
The theoretical value of the blend lifetime is 0.04 ns, 
which is so short compared with the time of the first data 
pOimcy OL they cascades ($33ns) >; > that® thew erteet of" the blend 
is probabably negligible. 

The lifetimes obtained from the ANDC analysis are shown 
in Table 4.4. The overall average of these values gives a 
lifetime of 0.32+.03 ns which is considerably lower than the 


values obtained from HOMER (0.416+.028) and TROY (.43+4.07). 
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Table 4.4 Results of ANDC analysis on I VII. 


PCC PANELING A PANELING B PANEGNG <C 
# T Come Ve = See ys se oh eG 
111 oO) «2 OWI 2a lt 22400. 020 0.75 wont ee 100 1.5% 
12:4 Mee a OVOP de ER. SAG 2024 3.5 oo Gee Oo 2 I 
211 sos .OtS As 7 $35 Pie lee 2343; 20048) 15:8 
Ze aoe Ol 2s 6 a SOOO e Sk eS DO Ues0 1 2a OX 


A Ve these tie 16 21526334 
PaNEUsg un .oybOr wt, 20:25) .30) 35 

B pe eo so lon e256 129.933 
PANG bopel “4 0°Oe 02) 16 20' 25° 32 35 

€ lee eee oe Vlri to 1 O21 Oe tae oi 5 Os he ois 
Prince |e cotene POUIZAIHE AG. 18020 22724426. 28: 30° 32 34. 35 


R--rejected due to large constant and/or chi-squared 

X--rejected due to large constant and negative transition 
probability 

*--initial average using best fit of each permutation 


* average is 0.32+0.02 (6.3%) 


An attempt was made to determine whether the blend in the 
464A cascade was having a significant effect on the lifetime 
by removing a single exponential term from the HOMER fit. 
The effect of removing the shortest lifetime term (0.24ns) 
Was to give essentially the same lifetime values but with 
higher reduced chi-squared values for some of the 
permutations (Table 4.5). The lifetime obtained from the 
weighted average of the values in Table 4.5 is 0.34+0.01 ns. 
This value is little different from the previous ANDC value 
Since the entire primary decay curve occurs at t>0.26 ns, 
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Table 4.5 "Blend corrected" ANDC analysis of I VII. 


PEC PANELING A PANELING C 
i Mion) oe T oS AADE 
tt or 2 Ore 0230) Soe O10" 1.s 
124 294000965106 Pye oh S101 oh vam 
210A OOS 6 VEE  te19 wee OOL Vent 
221 Foe ow Ot) 16.6.2 ast IRSOGO TEs 
average lifetime = 0.34+0.01 ns 


Most, OF the’ primary decay curve. Fig. 4.7 shows a typical 


example of the ANDC analysis for this data. 


Fig. 4.7 Sample analysis(M-D) for I VII: variables. 
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Removal of either the middle or largest lifetime 
component of the cascade resulted in negative transition 


probabilities and large constants for the ANDC analyses. 
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The theoretical lifetimes of the two direct cascade 
levels are about 0.07 ns which, when compared with the 
relatively long lived cascade decays (the decay curves were 
measured out to 4.4 ns), indicate that significant cascade 
effects should be seen in the primary decay. Thus one would 
expect a Significant difference between the ANDC analysis 
and the multi-exponential fit. 

Combination of both the "corrected" and "uncorrected" 
ANDewGesults yields, ja, laftetime of 0.3340.02 ns. The 
theoretical value! of the lifetime, 0.285 +<20%, agrees with 


the ANDC, result. 


The Se VI ion has the same number of valence electrons 
as I Vil which was examined in the previous section. Thus 
the decay scheme shown in Fig. 4.8 is the same as that of 
Fig. 4.6, except that the principal quantum numbers are 
lowered by one. However, in this case the cascade 
transition from the 2S levels is much weaker than that from 
the 2D levels. Furthermore, the data obtained for the 5s2S- 
Yp2P transitions indicate a very short lifetime (<.15ns). 
Thus only the one direct cascade from the -D° levels meed be 
included in the ANDC analysis. 

These data were taken with ANDC analysis in mind, so 


en ewe wwe a ee ee = 


1 Obtained by Kwok-tsang Cheng (University of Chicago) using 
Dirac-Hartree-Fock calculations. 
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Table 4.6 Characteristics of Se VI Data. 


PEAK ACCUMULATION 
(Ay NOs. SCOUNTS TIME (SEC) 
Primary 886 2 3000 60 
Cascade 588 Z 5000 20 
Primary 844 2 3000 30 (#1) /60 (#2) 
Cascade 606 2 4000 PAS) 


-last data point taken corresponds to t~11ns 
-50 points in each decay curve 
-data taken using multi-channel analysis technique 


Beam energy 1 MeV Beam velocity 1.55=mm/ns 


Fig. 4.8Selenium VI decay scheme - 4p*P levels 
(numbers in {} are approximate measured intensities) 
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the distance coordinates of the primaries and cascades are 
matched. However, problems with the foil stepping mechanism 
affected the distance measurements in some cases, So HOMER 
fits had to be used for the 886#1 and 844#5 cascades to 


compensate for this error. Table 4.7 shows the results of 
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the ANDC analyses using the (2-D), (2-A) 


and using three different partitionings. 


Table 4.7 Results of ANDC analysis on 
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partitioning C is considerably different than A or B in that 
there is a large degree of overlap between the adjacent 
panels and 13 points have been dropped from the tail of the 
decay curve. This was done due to the poor statistics of 
the last 13 points. The lifetimes obtained from 
Partitioning C are not significantly different from those of 
A or B, but there is a noticeable improvement in the reduced 
chi-squared of some of the permutations. This indicates 
that the choice of paneling does not affect the lifetime, 
but may affect the statistics of the analysis. Panes 4.9 
and 4.10 show a typical example for this analysis. 

The difficulty in analyzing the permutations with 886#1 
may be due to inadequate multi-exponential fits or simply 
pOCE cata Sets. “A\ Correction (of the beam velocity to 1.576 
mm/ns +2% lowers the lifetimes so that the average values 
are 0.35740.029 (844A) and 0.401+0.029 ns (886A). The ratio 
of the two lifetimes is 0.89+411%. Theoretically this ratio 
should egual the cube of the ratio of the corresponding 
wavelengths, which is 0.86. The resuits of TROY multi- 
exponential analysis yields the same lifetime for both 


Vevels: (055240201) ns, 


This chapter demonstrated the practical application of 
ANDC analysis to actual experimental data. An important 
consideration when performing an experiment is the 


integpation length.  Z£OANDC aseto bezapplicd.to thetdata, 
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Fig. 4.9 Sample analysis (2-A) for Se VI: graph. 
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Figs 4.10 Sample analysis(2-D) for Se VI: variables. 
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all decay curves that are to be used in a given analysis 


should either have the same integration length or the 


integration length should be so small as to have a 


negligible effect on all the decay curves, 


The lifetime of the 3p3D level in 


II was determined 


to be 6.994+0.34 ns (including the uncertainty in the 


velocity by Using ANDC analysis sis this analysis gave no 


reason to suspect any problems such as undetected blends or 


missing direct cascades and the result is also in excellent 


agreement with the conventional multi-exponential analysis, 
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one can be confident in the lifetime obtained. 

In the case of I VII, the lifetime of the ePoe sy oerevel 
was determined to be 0.3340.33 ns by ANDC analysis. 
However, even though the ANDC analysis is self-consistent, 
there 1s some doubt of the validity of this result due to 
the blend present in the direct cascade at 464A. This 
result is also significantly lower than that given by multi- 
exponential analysis, although it agrees with the 
theoretical value. Analysis of further data, with the blend 
at least partially resolved, is necessary to obtain a 
conclusive value for the lifetime of this level. 

The lifetimes of the 4p2P3/9 and 4p2P1/2 levels of Se 
VI were determined by ANDC analysis to be 0.3574+0.029 ns and 
0. 401+0.029 ns respectively. The analyses were such that 
confidence can be placed in this first value. The latter 
value is not as certain due to the problems of in one of the 
two primary runs analyzed (one should be wary of drawing 
conclusions from only one "good" primary). However, since 
the lifetimes have the correct theoretical ratio, the second 
lifetime is probably correct also. As¥in the previous case 
the ANDC lifetimes are much less than those obtained by 


conventional analysis. 
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CHAPTER V 
CONCLUSION 

The particular method of applying ANDC analysis 
developed in this work has shown itself to be a satisfactory 
means to determine lifetimes in decay situations involving 
Significant cascading. It has particular value where multi- 
exponential techniques fail to give satisfactory results, 
but is also useful as a check on the results of multi- 
exponential analysis. 

In the highly cascaded synthetic data situations of 
Chapter III, the multi-exponential method proved to be 
inadeguate, giving either misleading results (ie. actual 
lifetimes not within one standard deviation of the 
determined value) or results which were only correct by 
virtue of very large error estimates (+20%). The ANDC 
lifetimes on the other hand were accurate to well within 
reasonable error limits (<7h-in°%alivcases)is(One of ‘the 
synthetic data cases also showed that cascade effects will 
not necessarily cause the lifetime obtained by multi- 
exponential analysis to be too large, as had been previously 
accepted in the past. 

The analyses of the real data of Chapter IV indicated 
that there is a point of diminishing returns in including 
datatpoines from ithe ttdi beofdasdecaytcunve wherewthere sare 
fewer than about 100 counts per data point, since the 
statistical fluctuations will then exceed 10% of the signal, 


causing the error in the intensity differences, A, to be 
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very large. The weighting procedures in the fitting routine 
insure that these points will have a small effect on the 
lifetime, but their presence certainly does not improve the 
analysis. 

Both the synthetic and real data cases indicate the 
desirablity of using as many different primary runs as 
possible in the ANDC analysis. Since the primary data set 
dominates the determination of the lifetimes and error 
estimates, this practice produces a number of independent 
lifetime measurements which can then be averaged. Three 
primary data sets would appear to be the minimum number to 
be coilected for a given transition. In the case where one 
data set proves to be unsatisfactory in the analysis, this 
would still leave two data sets to provide two independent 
ANDC analyses whose results can be checked for consistency, 
and then averaged. Simi lariypirt would also be/wise ‘to 
obtain at least two data sets for each direct cascade. 

Before collecting the decay curve data, the decay 
scheme and the spectrum {at the energy at which the decay 
data is to be taken) should be studiedsecarefully to 
determine all the possible significant direct cascades which 
should be included in the ANDC analysis. At the same time, 
the possibility, of blends in any ofthese decays ‘should be 
considered. It is also necessary that all the decay curves 
to be used in a given ANDC analysis be taken at the same 
beam energy to insure that the level populations are the 


same for each decay curve (a basic assumption in ANDC 
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analysis). As indicated in Chapter IV, it is also necessary 
that the optical integration length along the beam be the 
same for all decay curves or that it be negligible for all 
decay curves in a particular ANDC analysis. 

The computer program developed allows the ANDC analysis 


to be performed several different ways (see Table 5.1). 


Table 5.1 Variations possible in the computer program used 
for ANDC analysis. 


a) FITTING PROCEDURES 


EQUATION VARIABLES NAME OF ANALYSIS 
P/A=1 +2 C/A Y=BR/Aph=C7A (2-D) 
A/P=1/t- &C/P Yabo (2-A) 
A=P/T-EE, C/t Z=A ,Y=P,X, =Cj (M-D) 


b) FORM OF DECAY CURVE DATA 


CURVE FORM DESCRIPTION 

PRIMARY -REQUIRES ORIGINAL HOMER NORMAL MODE 
FORMAT DATA SETS 

CASCADE -ORIGINAL HOMER FORMAT WITH NORMAL MODE 
MATCHING DISTANCE VARIABLES 

CASCADE -ORIGINAL HOMER FORMAT WITH LINEAR INTERPOLATION 
MIS-MATCHED DISTANCE VARIABLES (L-I) MODE 

CASCADE -SIX PARAMETER EXPONENTIAL MULTI-EXPONENTIAL 
FIT TO DECAY CURVE HOMER FIT (MHF) MODE 


*k*k note that all the cascades for a given ANDC analysis 
must be input in the same mode. 


c) Analysis can be done using any desired panel division 
of the decay curves. 
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Previous chapters have discussed the relative merits of each 
variation on the method. However, since all the variations 
are capable of producing a "valid" ANDC analysis, performing 
the same analysis using a few different variations is a 
useful check on the self-consistency of the results. 
Variations of the paneling does not usually affect the 
lifetime significantly but can give a better fit in the 
sense of lower error estimates and/or reduced chi-squared 
values. 

The results of an ANDC analysis may be accepted as 
valid subject to the satisfaction of the rejection criteria 
on various quantities resulting from the analysis which were 
discussed in sec. 3.4. Once an acceptable lifetime has been 
obtained for a particular level (following all the 
procedures outlined above), it is of interest to attempt a 
multi-exponetial fit using the determined primary lifetime 
and the cascade lifetimes (as determined by HOMER fits) as 
Tipeteal jestinates for ‘che line. (Che sec.ats.3)hp ipossrbiy. 
constraining one or more of the fitting parameters. 
Frequently, the lifetime determined by the multi-exponential 
fit will agree with the ANDC lifetime, indicating that the 
decay "situation wasijnoti/so sever ebyicdscaded.asiito render 
conventional procedures inadequate. This will also serve as 
a check on the error estimate of the lifetime. In the cases 
where the multi-exponential fit and ANDC lifetimes are 
inconsistent, "tie SANDG (liféetine *1s*toxbe wprefermed. The 


only reason an acceptable ANDC analysis would be incorrect 
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is when a Significant direct cascade had been omitted or an 
undetected blend had occurred in either the primary or one 
of the direct cascades. AS was discussed in sec. 3.4, there 
is no absolute guarantee that the application of the 
rejection criteria will always detect these possible 
problems. However, these situations are less likely than 
the possibility of some difficulty with the multi- 
exponential analysis (cf. sec. 3.3). 

Refinements in the application of the ANDC method can 
be made by further study of synthetic data situations. The 
relationship between blends and the determined values for 
the lifetimes, estimated errors and the rejection criteria 
can be studied by varying the type and amount of blend 
present. A less qualitative "curvature" rejection criterion 
might be developed by calculating the directional deviation 
of the fitted curv.e from the actual ANDC points. Also, such 
studies could place definite limits on what constitutes a 
"significant" direct cascade. The relationship between the 
number of counts (peak and minimum) and the quality of the 
resulting ANDC analysis could be easily studied using 
synthetic data situations. This would be of value in 
determining how many counts one needs and how far into the 
decay curve tail on must go in order to obtain an optimal 
data set for ANDC analysis. 

In the original paper describing the ANDC method 
(CBB 70), it was stated that the method was "very 


promising". However, since then, very little use has been 
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made of this technique due to an apparent lack of success in 
Tes practical application. ' The results of this work show 
that the ANDC method is able to fulfill its original 


promise. 
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APPENDIX I 


Modeling of Cascaded Decay Schemes 


The population of a particular level in a known decay 
scheme can be determined as a function of time if all the 
relevant initial populations and transition probabilities 
are known. This is done by solving the coupled differential 
equations governing the decays of each state, as was done 
for the three level system*in sec. 1.3. The equations have 
been solved and put in mnemonic form by L.J. Curtis (Cu 68) 
for a generalized decay scheme. The formula which follows 
refers to a decay scheme in which the energy levels are 
labelled from 1 to min order of increasing energy. The 


population of level n is then given by 


te m m-1 m Bi. 
eee Ne kO)ec 7 + pad thea!) ee ) ) {j>icn} + 
i=n+1 i=nt+l j=it+l 
CAT S15) 
m-2 m-1 m 
a ) ) Sete isin ee. Sew ine Peder, 


i=nt+l j=itl k=j+l 


where the brackets correspond to contributions from single 
stage, two stage, three stage etc. transitions into the 


state n, and 


Aatat a tnt 
{i-n} = N, (0)A, — a5 = 
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When modeling cascade decay schemes for the purpose 
of lifetime determination, it is more convenient to have 


the level populations expressed in the form 


NA (te), = )D e : 
ae 


This can easily be achieved by factoring eqn. Al.1 into 
j-l 
Fe * WL Aa+1 4 . fs f 
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k=n 
k#r 
Thus Lf one has a decay Sschemevas shown in) Firgq..Al. 1, 


one can then write for the intensities of the primary 


and direct cascades, 


m Seas a c 
Pit) = Peer and. C(t) %=) )we@%e , where 


P =A 


de . 
a oe and Cr= Az, Qrg - Le7siouldy be noted tthat. the 


index, Gd, “eters to a’ particularecascade path, So that. the 


same level may occur in more than one path. 
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APPENDIX II 

A description of the input parameters and their format 
are contained at the beginning of the listing for the main 
program. The two borrowed subroutines, REGRES and CORREG, 
and the simple plotting routine, APLOT, are not listed here. 

The program” output consists of a partial listing of the 
input parameters (not including data sets) and the results 
of the analysis, samples of which were given in Chapter IV. 

Phe results of (27D) or (2-A) analyses! consist of 
- the calculated ANDC variables(X(I) and Y(I)) with a plot 

of these variables (if desired), 


-the calculated variances and covariance (VX,VY,VXY) of the 
Variables, 


-the weight(W) calculated for each point, | 


-a measure of the deviation(D) of the fitted line from the 
calculated variables, 


=x (2) ,YO(1), the: fitted values forrthe varrables, 
-the center of "mass"(CENTROID) of the calculated variables, 


-the parameters of the fit (SLOPE and INTERCEPT), where the 
intercept is the lifetime (2-D) of its,inverse (2-A), 


=the standard deviation(S.D.) of the CENTROID, SLOPE and 
INTERCEPT, 


-the SQUARED RESIDUALS and # DEGREES OF FREEDOM, the ratio 
of which is the reduced chi-squared, 


-other output related to the details of the computation 
(MATRIX, LATENT ROOTS e€tc.). 


a es ee oe 


1 For a description of the mathematical details of this 
fitting procedure see CRRW 72. 
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The output of the (M-D) analysis consists of a check 
the validity of the assumption that has a significantly 
larger error than any of the areas, 

VAL=(V (A) /A2) / ees) aes and 

-the ANDC variables DELTA and AREAS , 

=the FITTED value for DELTA, 

-the LIFETIME and its estimated standard deviation, 


-the transition probability (PROBE?) modified by the 
efficiency for each direct cascade, 


=the value of the constant, K, and its estimated error, 
-the value of the reduced CHI-SQUARED. 

The foliowing pages contain a listing of the main 
program used for the ANDC analyis and the subroutine which 
does linear interpolation for mis-matched (distances) 


cascade data sets. 
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PROGRAM FOR ANDC ANALYSIS OF DECAY CURVES 


DEVICE 7 HOMER PARAMETERS TO RECONSTRUCT CASCADE(S) 
LINE N AMP,TL,ALIGN, ICOM (7F9.4,4A4) 
AMP=AMPLITUDES TL=LIFETIMES ALIGN= INITIAL DISTANCE OF 
CASCADE FIT 

DEVICE 8 DECAY DATA IN HOMER FORMAT (P'S THEN C!S) 
DEVICE 9 CONTROL DATA 

LINE 1 IND,ND,IFLAG,IM,1S,1IF,NC, NPERM, VEL 

(OT P27 31k, UT 1, Lae Po 1) 


EENEY 2 1VPI (2013) 

LINE 3. PF (2013) 

LINE 4/5/6 ID(1,1) (3012) 

LINE 5/6/7 ID(2,1) 

LINE 6/7/8 ID(3,1) IF APPLICABLE 

LINE 7/8/9 ID (4,1) IF APPLICABLE 
IND----HOMER DEPENDENT/INDEPENDENT/I NTERPOLATION-1/2/3 
ND----- DISPLAY PARAMETERS 


1--NO PLOT,SHORT STATISTICS 

2--PLOT,SHORT STATS. 

3--NO PLOT,LONG STATS. 

4--PLOT,LONG STATS. 
IFLAG--(2-D) /(M-D) /(2-A) --0/1/2 


IM----- NUMBER OF DECAY CURVES IN ONE ANALYSIS 
Is----- FIRST DATA POINT OF PRIMARY USED IN ANALYSIS 
IF----- LAST DATA POINT TO BE USED IN ANALYSIS 

WHEN IF=0: PI,PF ARE SPECIFIED BY USER 
NC----- NO. OF DECAYS OF EACH TYPE (IE. #P,#C1,#C2, #C3) 


NPERM--NO. OF PERMUTATIONS DESIRED MAX=40 
=0 =>AUTO PERMUTATION 


VEE-==-BEAM* VELOCITY > (HN/NS) 
el -S== INITIAL PANEL POSITIONS(MUST HAVE PI(1) SMALLEST) 
BR ras4= FINAL PANEL POSITIONS (LAST MUST BE LARGEST) 
ICOM---ONE COMMENT FOR EACH DECAY CURVE (MAX-_0) 
LD==-=F PERMUTATIONS: NOT NECESSARY IF AUTO-P=:8M ON 
CAUTTON: “4 OF POSSIBLE, PERNS HAY BE > MAK. (ALLOWED. 
TSEeLINE PRIMARY ID%S A wie.g 4 
2ND LINE CASCATE ID'S 2'B%2*3 
ETC. 


SUBPROGRAMS REQUIRED TO BE IN USERS FILE: 
REGRES CORREG APLOT INTER 


COMMENTS: INTERPOLATION SHOULD BE USED WITH ONLY ONE 
PRIMARY CAT «A +TIMNE* (UNLESS P* S@DISTANCES, MATCH) 


REAL*8 XX(20),YY (20) ,VX(20) , V¥ (20) , VXY (20) 

TeSUNTSUMK, SONY, SUMXY, SUMX2e,SUMY2),2,2X 21 ,X1,¥1 

REAL Y (80) ,D(80) ,TA(80) ,AMP (3) ,TL(3) ,7(80) ,¥FIT(20), 
1SDT (200) ,YA (400) ,DEL (200), SYA(400) ,SD(200) ,R(40) ,A(4), 
2SIGA(4) 

INTEGER ICOM (4,20) ,ILO(3),NPP(20) ,ID (4,40) , THOM (5, 3) 
1,PI(21),PF (20) -NPNP(5) ,NC(4) ,ICOML (4,4) ,ICOM1(4) 
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COMMON /IDISP/ND /AREA/YA, NPAN,NPNP 

EQUIVALENCE (XX(1),¥(1))-(YY¥(1),¥ (41)),(VX(1),D(1)), 
T(V¥ (1) -D(41)), (TL (1), TA(44)) , (ILO(1) ,TA(47)) 

1, (VXY (1) ,TA(1)), (AMP (1) ,TA (41) ) 
2, (NPP (1) net (Ws (NC (1) @NP), (NC (1) 7 NPNP (2))) 
3, (NC1,NC(2)),(NC2,NC(3)), (NC3,NC(4)) 

DATARTHUOM/ OHO’, fMER. ©, "DEPE*, (NDEN', *T ay 
Tae. (Uy NDE ee PEND RE AYENT ©! LINE’, 8AR 7 
egret ye POLA!,*TIONY/ 

CALL TRACER (-1) 

1000 READ(9,100,END=2000) IND, ND,IFLAG,IM,1S,I1F,NC,NPERM,VEL 

PEQ21)\=0 

IF(IF.NE.0)GO TO 73 


PANEL SIZES 


READ(9,106)PI,PF 
DO: ATARI =3)420 
IF(PI(J).EQ.0)GO TO 72 
71 CONTINUE 
72 NPAN=J-1 
G0. 10) 77 


AUTOMATIC PANELING 


73 P=1IS+4 
DO 74 J=1, 20 
PI (J) =P-4 
PF(J)=P 
P=P+5 
ER(P2GT.LF)GQ TO 75 
74 CONTINUE 
75 NPAN=J 
J=J+1 
DO 76 K=J, 20 
PTI (K) =0 
76 PF(K)=0 
77 WRITE(6, 150) (IHOM(K,IND) ,K=1,5),PI, 
NS=PF (NP AN) 
IDUP=NS-1-PI (1) 
NPAC=NP+4NC (2) +NC(3) +NC (4) 


AUTOMATIC PERMUTING OF PRIMARIES AND  aASCADES 


IF (NPERM.NE.0)GO TO 13 
K=0 

DO 60 I=1,NP 

DO 60 J=1,NC1 
DORe0sL=1, NC2 

DO 60 M=1,NC3 

K=K+1 

GO (1AK p=1 

ED (2QEK ped 
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75 


ED(37/K)=L 

ID (4,K) =M 

CONTINUE 

NPERM=K 

IF(NPERM.GT. 40) GO TO 2000 
GO TO 78 

DO’ 700 I=1,16M 

READ(9, 102) (ID (I,J) , J=1, NPERM) 
CONTINUE 

M=0 

NPNP (1) =0 


READ IN DECAY. CURVES 


1500 


: 


DO! 14 I=1, NPAC 

IF(I.LE.NP)GO TO 1500 

GO TO (50,1500, 3000) , IND 

READ (8, 103) (ICOM (J,1I) ,J=1,4),N,BGR 
READ(8,104) (D(J) ,¥(J) ,TA(d) ,J=1,N) 
po. “105219 

Y (J)=Y¥ (J) -TA(J) *BGR 


MATCH DISTANCE ENDPOINTS 


IF(I.GT.1)GO TO 5 
DS=D (IS) 

E20 

DOp2ak=uS, i 

L=L+1 

T (L)=(D(K) -DS)/VEL 
GO TO 6 

DGs8.K=1,N 

IF (D(K).EQ.DS)GO TO 4 
CONTINUE 

1S=K 

DFATSS EO 1)GOLTO 10 
E=0 

DOA KSLSAN 

L=L+1 

¥(Uys_Y (K) 

GO TO 10 


CALCULATE CASCADE AMPLITUDES FROM HOM: FIT 


50 


of 


READ(7, 105) AMP, TL,ALIGN, (ICOM(J,1) ,J=1,4) 
ALIGN= (ALIGN-DS) /VEL 

DOWS6 K=1,3 

TLT=TL (K) 

TLOKKP=1 

IF(TL(K)~EQ.0.0)GO TO 57 

AMP (K) =AMP (K) *EXP(ALIGN/TLT) 

GO TO 56 

I LO (K) =3 
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76 


CONTINUE 
DO 55 K=1,NS 

B=0.0 

DO 54 J=1,3 

IBB=ILO (J) 

GOTO (51,54 353) ,1BB 
TLT=TL (J) 

E=-T(K) /TLT 

IF (E.GT.-35.)GO TO 52 
ILO (J) =2 

GO TO 54 
B=AMP(J) * EXP (E) +B 

GO TO 54 

B=B+AMP (J) 

CONTINUE 

Y (K) =B 

MBETE(6 5452) (CON (K,I) ,K=1,4) , AMP, u.,ALIGN 
GOL TO# 10 

CALL INTER (N,DS,D,Y, ICOM1) 
pors00ttJ=1,-4 

ICOM (J, 1) =ICOM1 (J) 


c CALCULATE AREAS AND DELTAS 


10 


222 


un 


NL=N-NS 
WRITE (6,151) (ICOM(K,1) ,K=1,4) ,1S,. -,NPAN,VEL 
T=0 

DO 12 J=1,NPAN 

K1=PI (J) -PI(1)+1 

K2=PF(J) -PI(1)+1 

K2M1=K2-1 

AA=0.0 

SA=020 

DT1=7T(K1#1) -T(K1) 

DO th T3k1,/ 42" 

DT=T(L+1) -T(L) 

AA=AA+ (Y (L)+¥(L+1)) *DT 
DPESEOUK1)GO TO 11 
SA=SA+Y (L) *(DT+DTO) **2 

DTO=DT 

SA= (SA+DT1*¥DT1*Y (K1) +¥ (K2) *DT*DT) /4. 
M=M+1 

YA(M) =AA/2.0 

HFAPSGT.NP)GO TO 12 

Q=0. 

SDP=0. 

IDT=1 

TD=T (K1) 

KC=K1-3 


Cc DOSLEASTVSQUARESSFIT FOR DELTA 


Pag 


SUM=0. 
SUMX=0. 
SUMY=0. 
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SUMXY=0. 
SUMX2=0. 
SUMY2=0. 
1G0=0 
PRK. UP. SAND IDP.EO. 1) IGO=2 
IF(K2.GT.IDUP.AND.IDT. EQe-1) IGO=-2 
221 Doe250- wrr=i1,'5 
L=KC+II+IGO 
XT=T(L) 
Tv (1) 
YI=DLOG (Z) 
ZY=Z¥*YI 
TX=ZEGI 
SUM=SUM+Z 
SUMX=SUMX+ 2X 
SUMY=SUMY+ ZY 
SUMX2=SUMX 2+ ZX¥*XI 
SUMXY=SUMXY+ ZY*XI 
SUMY2=SUMY2+ZY*YI 
250 CONTINUE 
DIV=SUM*SUMX 2-SUMX*SUMX 
CA=(SUMX2* SUMY-SUMX¥*SUMXY) /DIV 
Z=SUMX Y* SUM- SUMX*SUMY 
CB=Z/DIV 
JI=II+1 
R (JJ) =DABS (Z/DSQRT (DIV* (SUM*SUMY2-SUMY*SUMY)) ) 
CALCULATE DELTA 
251 P=EXP (TD*CB+CA) 
SDP=SDP+P/ (R (JJ) **2) 
Q=Q+IDT*P 
LF (IDT.LT.0)GO TO 260 
TD=T (K2) 
IDT=-IDT 
KC=K2-3 
GO TO 211 
260 CONTINUE 


DEL (M) =Q 

SD (M)=SDP/4. 

SDT (M) =SQRT(SD(M)) 

R (J) =(R (Jd) +R(JJ-1)) /2- 
12 SYA(M)=SA 

IF(I.GT.NP)GO TO 14 

I=507/2 

WRITE(6, 154) (R(L) ,L=1,JJ) 
14 CONTINUE 


IF(IFLAG.NE.1)GO TO 17 
{M-D) ANDC ANALYSIS 


DO 18 I=1,NPERM 
WRITE(6, 157) 


77 


4 at a) 14), ) : a : 7 : 
6S SOOT CF atid «/ DU ai r i, * ts a 
me, ar ay 7 i 
r 


J 7 
= it ; cf me eh ; Qt aA at. ares a 
PPh fats Ort OC 
a Omg t Shed a 
, my rz a 
7 6h, as5 nt 
 (entnaera aT 
ie : af 


‘> wae ees 
TREES _ 
AUEROE Lie 
wo r OLE. Te 
just aeeeane A 

ax eee NISASEABS 1 
hos Oa os. Ga 
Vai ¢ 4 wae ten ue P, : 
as Oe 
} ty. Mild pri eee oa 
TaN TRH A EOF he RE RN 
:< Reh eg ee enter We wilt 
ES ea 

‘ie . “puet: 
((r hileringe -S2 NGG ae Beyer td) ‘AORt\ Sh eis tay 
ae ERS STAG 


» 


a Ry A 
(i Teal Wy yay et. 

| ee tS 
he bo Gt Om ‘0. Rd ROT $e Sd 
| maa Sass 


Poa ne 
‘i 


| oR 4 a a icky Pee ay Ei - 
| gee ieee 


‘J ¥ 8 xs or oe hae 


+ 
an 
re 
eo 


7 


rn 
= 


78 


L=NPAN* (ID(1,I) -1)+1 
CALL REGRES(A,DEL(L) ,SDT(L) ,NPAN,IM,ID(1,1),1,YFIT, 
1A0,A,SIGAO,SIGA,R,RMUL,CHISOQOR, FTEST) 
IN=0 
DO 23 M=1,I1M 
NN=ID(M,1) +IN 
IN=IN+#NC (M) 
DOee4aJ=1, 4 
24 ICOML(J,M)=ICOM(J,NN) 
WRITE(6, 160) (ICOM(J,NN) ,J=1,4) 
23 CONTINUE 
160 FORMAT(' *,4A4) 
WRITE (6, 162) 
162 FORMAT(*-","* VALIDITY" ,2X,*FITTED', <, 
1*S.DEV ', 4(4X,"AREA', 4X," VARIANCE) 
BOD=0.0 
DO 20 J=1,NPAN 
AOD=0.0 
DOp2ds K=44 IM 
M=NUM(Jd,K,1ID(1,1)) 
Y (K) =YA (M) 
D (K)=SYA(M) 
21 AOD=AOD+D (Kk) / (Y (K) **2) 
K=L+J3-1 
AOD= (SDT (K) /DEL (K) ) **2/AOD 
WRITE (6,161) AOD, YfIT(J) ,DEL(K) ,SDT(K) , (Y(K) ,D(K), 
qK=125M) 
161 FORMAT (1X, F6.1,1X,3(2X,F6.0) »8(2X, F8.1)) 
20 CONTINUE 
TAU=1. /A (1) 
DTAU=SIGA(1) ¥*TAU*TAU 
DO 1600 J=2,1M 
1600 A(J)=-A(J) 
WRITE(6, 164) TAU, DTAU, (A(J) ,SIGA(J) ,J=2,14) 
164 FORMAT (///! LIFETIME=",F7.3,1X,F7.3,16X,' PROBE? ',4X, 
‘T*ERRORS ',3(/40X,F8.4,3X,F8.4)) 
WRITE (6,165) AO,SIGAO,CHISOR 
165 FORMAT (‘OCONSTANT= ',F9.4,3X,F9.4,5X,'CHISQR=! ,F8. 3) 
WRITE(1, 183) ((ICOML(J,K) -J=1,4) ,K=1,1M) 
WRITE(1, 180) TAU, CHISOR,AO, (A(J) ,J=2,14M) 
WRITE(1, 182) DTAU,SIGAO, (SIGA (J) ,J=2, IM) 
180 FORMAT (16X,F7.3, 5X, F6e 3, 5X ,F904,5X,3 (F8. 4, 1X)) 
182 FORMAT (16X,F7.3,16X¢F904,5X,3(PS804,1X)) 
183 FORMAT(*® *',20A4) 
18 CONTINUE 
GO TO 1000 


‘DELTA % 2X, 


(2-D) ANDC ANALYSIS 


17 DO 16 I=1,NPERM 
L=NPAN*(ID(1,I) -1) 
M=NPAN* (ID (2,1) +NP-1) 
DO 15 J=1,NPAN 
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TP=i+I 
IC=MtJ 
C 
C ALTERNATE (2-D) METHOD 
E 
IF (IFLAG.EQ.0)GO TO 8000 
O=YA (IP) 
D2=Q*Q 
SDP=SYA(IP) 
AOD=DEL (IP)/Q 
VY (J) = (SD (IP) #+AOD*¥AOD*SDP) /D2 
GO TO 8001 


8000 SDP=SD (IP) 
Q=DEL (IP) 
D2=Q*Q 
AOD=YA(IP) /Q 
VY (J) = (SYA (IP) +AOD*AOD*SDP) /D2 
8001 BOD=YA(IC)/0 
VX (J) = (SYA (IC) +BOD*BOD*SDP) /D2 
VXY (J) =AOD*BOD*SDP/D2 
YY (J)=AOD 
15° XX (J) =BOD 
T1=1D(1,1) 
T2=1ID (2,1) +NP 
WRITE(6, 153) (ICOM(K,1T1) ,K=1,4) , (ICOM (K,T2) ,K=1,4) , 
11S,1IF,NPAN 
GO TO (601,602,601,602) ,ND 
602 CALL APLOT (NPAN, XX,YY) 
WRITE(6, 157) 
601 CALL CORREG(NPAN,XX,YY,VX,VY,VXY) 
16 CONTINUE 
GO TO "1000 
100 FORMAT (411,213, 1X,411,13,F 6.41) 
102 FORMAT (4012) 
103 FORMAT (4A4,54X,12/16X, F6. 3) 
104 FORMAT (5 (F3. 1,F6.0,F4.1,1X)) 
105 FORMAT (7F9.3,4A4) 
106 FORMAT (2113/2013) 
150 FORMAT(*1",10X,5A4//10X,' PTS. IN PANELS ARE ° 
1/2113/2013//) 
15 TEFORMAn CTOX SYAW PSTS=" 13,2, "ILF=",13, 2%," NPAN=",13, 2%, 
1*VEL=',F5. 3) 
152 FORMAT (10X,4A4," AMPLITUDES',3F10.2,' LIFETIMES', 
203988. 3) 3X, “ALIGN=", F552) 
153 RORMAT(' 1,404," AND ",4AK IM TS=',173,2%,'1F=",13,2X, 


2'*NPAN=',13,13//) 
154 FORMAT(* LIN CORR. COEFF.'/1X,20('*",F5.2)) 
157 FORMAT('1®*) 
2000 STOP 
END 
FUNCTION FCTN(Z,1,d,NX) 
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DIMENSION NPNP(1),YA(400),NX(1) 
COMMON /AREA/YA,NPAN, NPNP 
M=NUM (I,J, NX) 

FCTN=YA (M) 

RETURN 

END 

FUNCTION NUM (I,J,N) 

DIMENSION N(1),YA(400) ,NPNP(5) 
COMMON /AREA/YA,NPAN, NPNP 

M=0 

K=n {5} 

DO* qhLs4A5 

M=M+NPNP (L) 

NUM=NPAN* (K#M-1) +1 

RETURN 

END 

SUBROUTINE INTER (NP,DS,X,YC, ICOM1) 


THIS PROGRAM READS THE CASCADE DATA SET AND 
CALCULATES THE AMPLITUDES 

Ate THSESDISTANCE COORDINATES GIVEN BY THE PRIMARY 
(LINEAR INTERPOLATION) 


REAL TA(80),YC(80) ,X(80) 
INTEGER ICOM1(4) ,ICOM2(4) , IX (80) , IXC (80) ,1YC(80), 


1IYS (80) 


READ(8, 102) I1COM1,NC,ICOM2, BGR 
READ(8, 103) (IXC(J) , LYS (J) ,TA(J) ,J=1,NC) 
POS2005=1, NC 


20 IYS{J)=1YS (J) -BGR*TA(J) 


PIND- PIRSTVUPSEMARY DISTANCE (POINT 
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1 


DON 24 fre 49 NP 

IX (I) =X (I) *10.0+.001 
IS=DS*10.0 + .001 
DOBIAIT=47NP 
PEGI 1). B0.1S).G0) 10 (2 
CONTINUE 


INITIALIZE INDICES 


SCAN CASCADE FOR DISTANCE POINT «GE. PRIMARY POINT 


SDS itA 


iC=1xC (1) 
LE (lec). NorGOwtOn 10 
TF (LC-1S5), 00,4 
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AAPIEG AC. eee S7kae 


wD 


5 ae 
rea 


Cra 
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OC} 


WHEN DISTANCE ENDPOINTS EQUAL 


30 J=J+1 
IYC (J) =1YS (I) 
GO TO 6 


HHEN NOT EQUAL DO LINEAR INTERPOLATION 


HY 2=175 (2) 
Y 1=IYS (I-1) 
KO= Agel) 
<= Tes (de 1) 
5 Q=I15S 
S= (¥2-Y1)/ (X2-X1) 
B=¥ 1-S*X1 
¥ 1=S*Q+B 
J=J+1 
IMC) =e 0 


ChEC Ka Dra NEXT P RIMARYs DIST. LSeWLIAINCURRENT 
CASCADE POINTS: ITF S50; 
USE SAME POINTS FOR INTERPOLATION. 


6 K=K+1 
IF(K.GT.NP)GO TO 10 
ms- Eh) 
IF (IC-IS) 3,30,5 


PINTS a ORF 


10 NC=J 
Dowell =n 
14. YC(i)=1YC(1I) 
URETE (6, 10M) (ox (2)., BYC(1) a =A ele) 
101 FORMAT (5(1X,13,16)) 
102 FORMAT (4A4,54X,12/4A4, F6.3) 
103 FORMAT (5(13,16,F4. 1,1X)) 
2000 RETURN 
END 
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APPENDIX ETL 


Invariance of ANDC Analysis under 
Time Integration 


Measurement of experimental decay curves requires 
collection of the emitted light.over,a.finite,.section of 
the beam. This effective time integration of the decay 
curve can cause a Significant deviation in the shape of 
the measured decay curve from the theoretical shape, depend- 
ing on the relative sizes of the averaging length and the 
lifetime components of the decay curve. Due to the nature 
of ANDC analysis, although the variables involved (ie. areas 
and differences in intensity) are affected by this integra- 
tion, the parameters obtained from the analysis (ie. slope 
ano intercept). ane, ain principle, not affected by the value 
of the integration length, as will be shown in the following 
paragraphs. 

If the “true" decay curve is given by 


T,(t) = ew eral aree ee : (Koay) 


then the measured curve when the integration length is 6 


will be 
t+6/2 


LG) = A cette ete! ‘ 
t-6/2 


Substitution of A3.1 into this then yields 
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tile. Jt SLOLet 10720 anadr where 
the common factor, 6, has been incleded in the detection 
erervecilency,( Ee for (chet transition. 

First, consider a four level decay scheme with only one 


direct cascade, as _ shown in Fig. A3.1.  ANDC theory states 


Fig. A3.1 Four level decay scheme. 
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condition can be re-written, dropping the (I,F), as 
ee (DeiGs) ese NED tae Cn) (A3.2) 


Expanding the quantities in this last expression into 


a sum of exponential terms yields 


Pp = Eipibi tet tre PB yoy + Be (epacn Bese iar 


fy ee ers ae ) faz} 


= E19 (Py2y/01 + P5Z5/o2 + P.2,/a3} ; 
C= By {Cn Zo/ca ata C32,/a3} 1 


I E49tP 121 J Bae + PZ} , mes 
a 

p= E491 P48 424/01 + P5S525/d2 + P48,2,/a3} ; 

Ck= EB, 1{C58525/o2 + C,842,/a3} : 


A'= Ej9{P,812, + Po8525 + P38325} - 

The P. and Cc. may be expressed in terms of lifetimes, 
transition probabilities and initial level populations, as 
was discussed in Appendix 1. In order to remove the consid-= 
eration of efficiency, one can write € = QEj9/Eo1 , where 
q = Ajot . The condition for the validity of ANDC analysis 


(eqn.A3.2) can now be written as 


{P,S,2,t+P.S525tP48,2,}{P12,/o1 +Z, (P4-qC5) /o2 + 


+2, (P-qC,) /as} = 


={P,2Z,+P2Z5+P32,1{P,S, 2) /o1 +85 Zo (Py-GCy) /o2 +823 (P,-GC,) /as} 


PPP (2 hed con Opn 3) 
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Note that, the efficiency factors, E and E, have cancelled 


10 Ace 


out. 
This equation can be most easily verified in parts by 
equating the factors of each =: 
So: P5Z5{P,2,/o2 + Z(P,-qC3) }= (P4-qC5) Z5{P2,+P 323} /a2 
3 


It can be easily shown(see below) that in general 


Pasdc = PO p/O1 Z (A3.4) 


Therefore, the previous 


equation becomes 


Po{P,2,/o1 + PZ} 
1 
Om ih 


The other two factors give 
to be invariant under time 
In order to prove its 


first must be shown that 


eats: 
Ee qc, 


d 
(hseis 1/5) 


& P5{P 12, + P4Z,}/a1 


= 1, 


the same result, showing ANDC 
inceqration, i1n thewrour level case. 


validity for any decay scheme, it 
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Using the results of Appendix I, in general one can write 
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Therefore, 


Aa nee 
Cae = Agi/ Ajo Cai-a) a Aio é cuales 
AyouCy = On/ 01  7eProving, eqn. WAS 4a. 


1 - qc. /P_. =]1- aoe 


In the case of the general decay scheme, eqn. A3.2 


becomes (note that C,=0), 


1 


{ }P,S,2.}{}2, (P,;-ac,) } = {}P,2,}{)8,2, (P,-aC,) /a, : 
ab 1 STC 1 1 
Oy 
Substitution of eqn. A3.4 yields, 


{)P,S,2,}{)2.P,/oz} = ()P,2,}{)8,2,P,/o,} . 
i i i i 
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